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transits 
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ABSTRACT 

We investigate gravitational lensing effects of an extrasolar planet transiting its host 
star. We focus on the 'rising spikes' of the light curve just before and after the transit, 
which is a peculiar feature of the gravitational lensing, and find that it could be a 
novel observable for determining physical parameters. Detectability of such an effect 
is also discussed. 
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1 INTRODUCTION 



Observations of decrease of starlight by a transit of an ex- 
trasolar planet is one of the important method not only for 
finding the planet itself but for determining physical param - 
eters of the transiting planet (|See e.g. Cassen et al.ll2006h . 
Roughly speaking, there are two observables in the transit 
light curve: Duration and depth. The latter gives the planet 
radius if the star radius is known, while the former leads to 
the planet orbital radius if one neglects the size and mass of 
the planet for the transit duration and the orbital velocity 
of the planet, respectively. 

There is often observations of the orbital velocity and 
the period of the transit using the Do pplcr ef fects in add i- 
tion to the transit observation (|See e.g. Cassen et al.|[200^ . 
Hence the planet mass can be determined uniquely, since 
the inclination is known to be near edge-on from the transit 

To date, hundreds of the extrasolar planets are observed 
by various methods mostly by Doppler method and tran- 
sitfl but those with a large orbital radius more than ~10 
au are limited by these methods. This is because one could 
see the transit practically only once, at most, in one's life 
time, and such planets are also difficult to detect with the 
Doppler method. 

As the distance between the star and the planet be- 
comes larger, the mass of the planet could in principle af- 
fect the transit light curve due to microlensing effect. Mi- 
crolensing with occultation has been cons idered generally or 
in other context { Macdcr 1973; Bromlevl[l 996: Marsh 2001; 
Bozza et al.ll2002l: lAgoUOOa : .Beskin fc Tuntsov.,2002. : Agoli 
20031 : ISahu fc GillilandlbOoS iLee et al.ll2009l ). It makes the 
depression shallower, and steepens the slope of the ingress 
and egress of the dip (,Marsh .20011) . However, the former 



can be explained by a smaller planet radius, while the latter 
could be due to a faster transiting velocity. 

A peculiarity of the microlensing effect on the transit 
light curve would be 'rising spikes' just before and after the 
transitUrhis is what we focus on in the present article, and 
we show that the height of the spike has different param- 
eter dependence from other observables such as the depth 
of the light curve. In particular, it could be measurable in 
high precision observations such as by Kepler or other future 
facilities. 

The structure of the article is as follows. In the next 
section, we consider both microlensing and occultation si- 
multaneously to obtain the magnification of the starlight. 
We see typical observables from the transit light curve, and 
define a novel observable as the 'rising spike' in the lensed 
case in Sec. 3, so that we consider three observables: The 
depth, the duration, and the rising spike. In Sec. 4, we show 
that the section of isosurfaces of these three observables will 
determine the physical parameters of the planet: The mass, 
the radius, and the orbital radius. In Sec. 5, we include the 
limb darkening as an example effect which may obscure the 
rising spike. We conclude in Sec. 6. Appendices are devoted 
to the other sections of the isosurfaces of three observables, 
and the other CEises for lensing and occultation in the aligned 
situation. 



2 MICROLENSING AND OCCULTATION 

The lens equation for a point mass is written as 



m 



(1) 



* E-mail: kasuya@kanagawa-u.ac.jp 
^ Sec http : //exoplanet . eu/. 



^ This feature was mentioned f or the case which occultation and 
microlensing effects are equal in ISahu fc~Gil liland 1,2003! ') , and in- 
vestigated for large source case in I Agoll 120031 ) . 
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Figure 1. Schematic view of the geometry of the gravitational 
lensing. 



where 8 is an angular separation between the lens and the 
image of the source, while /3 is that between the lens and 
the source without lensing (See Fig. 1). Be is the Einstein 
angle, 



AGMl Di 



DosDoL ' 



(2) 



which represents an effective size of the lens. Here Dos and 
DoL are the distance from the observer to the source and 
the lens, respectively. Dls is the distance between the lens 
and the source. Ml is the lens mass, G is the gravitational 
constant, and c is the speed of light. Solving Eq.([l|, we ob- 
tain the solutions 



0± 



(3) 



where 9+(-) is the angle of the image whose light comes 
from the far (near) side of the lens, which we will call it 
+ (— ) path hereafter. 

The magnification of the light of a point on the source 
star through the ± path is given by 



A+ = 



dl3l3d4> 



2/3 



■2 + 4^1 



(4) 



since the surface brightness is conserved. In the case of mi- 
crolensing, one cannot discriminate these two magnified im- 
ages separately, so the observed magnification should be the 
addition of the two as 

P'^ + 29% 



A{l3) = A+ + \A-\ = 



/3^/32+44 



(5) 



Since the source is usually an extended object, integrating 
over the source, we can obtain the total magnification as 



A- 



area with lens 
area without lens 



/ A{l3)l3dpd<j} 
J l3dl3d(P 



(6) 



So far we have considered only the point mass lens. 
Because we are interested in the observation of the transit 
of a planet in front of the source star, we need to investigate 
the finite size effect of the lens planet. There are three cases: 
(1) Both -I- and — paths of the light from a certain point on 
the source star are occulted, (2) — path is occulted, or (3) 
neither path is occulted. Therefore, when integrating over 



the source star in the numerator of Eq.(|6]), we must use 
A{I3) itself for case (3), A+(P) for case (2) instead of A(/3), 
or do not integrate over the region where case (1) holds. 

In the simple configuration that the centres of the lens 
planet and the source star are aligned, one can calculate 
the magnification very easily0 The denominator of Eq.® is 
estimated as 



(7) 



as usual. Here 9$ ~ Rs/Dos, where Rs is the radius of 
the source star. In order to calculate the numerator in the 
transit case when the occultation occurs, we must specify the 
condition to hold. The occultation takes place if I3{9l) > 0, 
where Ol ~ Rl/Dql with Rl being the radius of the lens 
planet (See case I in Fig. lAlfl . It leads to Ol > 9e, which 
implies that the size of the planet is bigger than that of 
the lens. Therefore, there is full occultation in the region 
13 < I3{9l), while - path is shielded for /?(6iz,) < P < 9s. 
Thus the numerator of Eq.([6} is obtained as 



A+{l3)/3dl3d(l) 



,8(eL) 

TT 
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fs~i3l + esy/9i + 4e% 



where Pl = P{6l)- Therefore, we have 



A- 



1 + 401 



Pl 



Pl + Ae% 



(8) 



(9) 



where the tilde denotes the variable normalised with respect 
to Os such that Pl = Pl/9s and 9e ~ 9E/ds- Notice that 
this magnification (or, more precisely, the diminution) cor- 
responds to the depth at the centre of the transiting light 
curve. For other situations, see Appendix A. 

In order to follow the whole process of the transit, we 
integrate numerically the numerator of Eq.® to estimate 
the magnification. We show some example of the light curve 
in Fig. [H Here we take Ml = lOMj, Rl = Rj, and a = 
200 au, where Mj and Rj are the Jupiter mass and radius, 
respectively. For parameters of the source star, we set, here 
and hereafter, Ms ~ Mq, Rs = R©, and Dos ~ 10 pc 
when calculating numerically. We should mention that Dos- 
dependence is very weak for Dqs ^ Dls- Also notice that, 
in this article, we restrict ourselves to the edge-on case in 
order to focus on the effect of lensing. 



3 OBSERVABLES FROM THE LIGHT CURVE 
AND PHYSICAL PARAMETERS 

3.1 Unlensed transit 

During the transit a light curve of a star drops down due to 
the occultation. We can see how much it drops down and 
how long it lasts. These are the depth and duration of the 
transit, respectively, and the only observables we can obtain 
from the light curve (in a simplest case). If we measure the 
depth at the center of the transit, it is estimated as 



More g eneral cases can be analytically estimated using elliptic 
integrals llAeolll2002h . 
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Figure 3. Schematic view of the hght curve of the lensed transit. 



Figure 2. Light curves of the transit. Red (solid) and green 
(dashed) lines represent the lensed and unlensed cases, respec- 
tively. The vertical axis denotes the magnification A, while the 
horizontal axis shows the time in day from the centre of the transit 
in the case of the circular orbit. We take Mi^ = lOMj, i?£ = Rj, 
and a = 200 au for the lens planet, and Ms = Mq, Rs = R©, 
and Dqs = 10 pc for the source star. 



Dos 



(10) 



A = ^ — ^OS 

where a = Dls is the orbital radius of the planet around 
the star. Once the parameters of the star [Rl and Dos) are 
known, it could be a function of Rl and a: A = A(_Ri,a). 
On the other hand, the duration of the transit is determined 
as 

_ 2{es + eL) _ 2{es + 0L) 



2{Rl + 



vl/Dol 
Dos 



Dos 



Rs 



G{Ms + Ml) ' 



(11) 



where the planet is assumed to revolve in a circular orbit 
around the mass center of the star and t he planet, and has 
orbital velocity vl = \J G(Ms + Ml) /a. Again, if the pa- 
rameters of the star [Ml, Rl, and Dos) are known, it be- 
comes a function of Ml, Rl, and a: t — t{Ml, RL,a)- If 
Ml could be much smaller than Ms to be neglected, we 
can derived physical parameters Rl and a from the transit 
observables A and t. Usually the transit observation is not 
the only information in our hands. Orbital velocity and/or 
orbital period of the planet is obtained from the Doppler ob- 
servation, and three physical parameters are all determined. 



3.2 Lensed transit 

Neither orbital velocity nor orbital period has been observed 
by Doppler method if the planet rotates around the star 
with relatively large orbital radius, typically farther than 
~10 au. In such a situation, the gravitational lensing effect 
becomes larger as the lens planet goes away from the star 
(See Eq.©). 

The gravitational lensing is known to affect the light 
curve of the transit in such a way that the depth becomes 
shallower and the slopes of the ingress and the egress be- 
comes steeper, but the former effect on the light curve can 
be mimicked by the smaller planet radius, while the latter 
may imply the faster orbital velocity of the planet. There- 
fore, one cannot tell if it is really affected by the gravita- 



tional lensing. In order to confirm that one must surely take 
into account the gravitational lensing effects, their peculiar 
feature should be observed in the light curve. 

Here we claim that the 'rising spike' is the smoking gun 
of the gravitational lensing effect. It is the rise in the light 
curve just before the ingress (or just after the egress) of the 
transit, shown as ^ in Fig. [3] There are several key points 
that the rising spike is a good observable for determining 
the physical parameters of the planet. First of all, it is a 
unique consequence of the gravitational lensing effect, and 
cannot be mimicked by others. The effect could be observed 
in high precision observation such as the Kepler, especially 
when the orbital radius is relatively large, where the orbital 
velocity cannot be observed. For parameter determination, 
it has different parame ter d ependence from the depth or 
duration of the transit (|Agolll2005 V Therefore, it is quite a 
unique observable. 

The amplitude of the rising spike can be calculated from 
the magnification of the light at the moment when the lensed 
shadow of the planet has just been tangent to the star as 
shown in Fig. [l] It is thus given by 



i + C = - 



A+Cp{p,^))pdpdi;, 



(12) 



where all the angular separations are normalized with re- 
spect to 9s such as a — a/9s, and 



A+{P) = 



with 



1 + 



?2 + 26*1 



13 



Pi +131 

/3l + 1 + pcosifj, 
psmifj. 



(13) 



(14) 

(15) 
(16) 



Since the image of the source star is not occulted 
through + path while fully occulted through — path, we 
can easily calculate ^ using the magnification of the -I- path 
for a point mass lens (|Witt fc Madll994 lAgolll2002l ) as 



47rV(/?L + 2)2(4^1 +/?i) 

x{(/3i + 2)2(4^1 +/3!)i5(fc) 
~ [PICPl + 2f + wICpI + 2Pl)] K{k) 
+4/3!(el + l)n(n,fc)}-i. 



(17) 
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Figure 4. Geometry of the planet and the star at the moment 
of the peak of the rising spike. 



where K(k), E{k), and n(n, k) are the complete elhptical in- 
tegrals of the first, second, and the third kinds, respectively, 
and 



n = 1 



0: 



k = 



4^1 



(18) 



Since 9e is a function of Ml and a, and Pl is a function of 
Rl, Ml, and a, the rising spike ^ is a function of Rl, Ml, 
and a: £^ — ^{Rl, Ml, a). 

On the other hand, the depth and duration of the transit 
are modified due to the gravitational lensing effect. From 
Eq.([9ll, the depth is now estimated as 



i + 1 + 491 + pL\Pl+4ei 



while the duration is given by 
_ 2(1 + /3l) 



.(19) 



(20) 



Therefore, these two observables are functions of Rl, Ml, 
and a: A = A{ML,RL,a) and r = T{ML,RL,a). Notice 
that they reduce to the unlensed case pO|l and ((TT}, re- 
spectively, for vanishing lens mass, A'Il ~^ 0, which implies 
6e ^ and Pl Ol- 



ISOSURFACES OF THE DEPTH, THE 
DURATION, AND THE RISING SPIKE 



Now we have the depth p9[) . the duration (|20p . and the ris- 
ing spike (|17p , which depend differently on the three physical 
parameters of the planet such as the mass Ml, the radius 
Rl, and the orbital radius a. If these three observables are 
measured, we can reversely use these equations to deter- 
mine three physical parameters. Although it could be done 
in principle, the inverse problem is practically difficult. How- 
ever, we could numerically obtain the isosurfaces of A, r, 
and ^ in three dimensions of Ml, Rl, and a. Intersection of 
these three surfaces must be a point, the solution that we 
need. 

Figure [5] shows sections of the isosurfaces of A, r, and 
^ for Ml = WMj, 3Mj, and Mj from the top panel to the 
bottom, respectively. Red, green, and blue lines represent 
the sections of isosurfaces of the rising spike ^, the depth 
A, and the duration r, respectively. Shaded region denotes 
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Figure 5. Sections of the isosurfaces of the rising spike ^ (red), 
the depth A (green), and the duration r in days (blue) for 
Ml/Mj = 10, 3, and 1, from the top to the bottle panels, re- 
spectively. Shaded region represents where the hollow disappears 
A < 0, which we do not consider. 



where A < 0, which we do not consider because we focus on 
the transit observation. Since their parameter dependences 
look quite different in the section with constant Ml, one 
could derive the physical parameters Ml, Rl, and a from 
the figure. For other sections, see App. B. 

As can be seen in Fig. [S] the rising spike is as large as 
~ 10"*, which can be measured in high p recision observa- 
tion such as the Kepler (|jenkins et al.|[201Cll ) , for Jupitar-like 
planets with the orbital radius of a few hundred au . On the 
other hand, ^ = 7 x 10"^ - 7 x lO"'' for a = 100 - 1000 au for 
those super-earth planets with M ~ lOM^ and R ~ 37?®. 
Therefore, it has no hope to see the rising spike in these 
cases. 

In order to observe the rising spike, it must last long 
enough in addition to the height larger than 10"'' . Actually, 
the duration of the transit is about 8 days transit, while the 
rising spike becomes half the maximum height in 4 hours for 
Ml = lOMj, Rl = Rj, and a = 200 au (See Figs. [2] or E]). 
Therefore, it could well be observed by Kepler with the long 
cadence mode, which has the time resolution of 30 minutes. 
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Figure 6. Limb darkened liglit curves. We take tlie limb darken- 
ing parameters as (ai, 02)= (0, 0), (0.18, 0.10), (0.40, 0.25), and 
(0.65, 0.35), respectively from the top to the bottom at the centre 
of the transit. Other parameters are set the same as in Fig. [2] In- 
figure is an enlargement of the rising spikes. Dashed line denotes 
yt = 1 + 10-4. 



5 LIMB DARKENING 

Realistically, the surface brightness of the source star is not 
uniform, but has position dependence such that the centre is 
brighter than the edge. This phenomenon is called the limb 
darkening of the star. Therefore, one must take this effect 
into account when any phys ical consequence is drawn from 
the light curve of the transit (lAgolll2002l . l2003l ). In particular, 
one may wonder if the rising spike would be round or totally 
disappear. In this section, we show that this is not the case. 

The effect of the limb darkening can be parametrize 
quadratically as 



f LD {cos fi) 



ai(l — cos fi) — 02(1 — cos fi) 



(21) 



where fi is an angle between the normal of the star sur- 
face and t he line of the sight, and ai and 02 are constants 
of O(O.l) l|Claret. Diaz-Cordoves fc Gimene jllQQSI ) . We can 
calculate the light curve of the transit including limb dark- 
ening effects from Eq.((Bl) with the limb darkening function 
H2ip being inserted into the integrand in both the numerator 
and the denominator. We show the transit light curve with 
the typical limb darkening in Fig. [B] Here we plot the case 
with (ai, 02)= (0.18, 0.10), (0.40, 0.25), and (0.65, 0.35) 
for Ml = WMj, Rl = 107?,/, and a = 200 au. One can 
see that the amplitude of the rising spike would diminish 
but still be large enough to be obse rved with th e precision 
of 10~* achieved by Kepler (Jenkins et al.|[2010l ) even for a 
large effect of limb darkening. 

A comment is in order. Gravitational lensing effects are 
also seen as a bump at the bottom of the light curve if the 
limb darkening is negligible. Although it is another peculiar 
feature of the lensing, the bump will totally disappear when 
the limb darkening effect is taken into account, contrary to 
the rising spike (Agol2003). We thus consider the rising spike 
as the unique apparent feature of the gravitational lensing 
effect on the light curve. 



6 CONCLUSION 

We have studied the gravitational lensing effects on the light 
curve of the transit which an extrasolar planet plays the role 
of the lens. The lensing effect becomes sizable only when the 
orbital radius of the planet becomes larger than ~ 10 au. In 
such situations, the orbital velocity and/or orbital period 
are hard to be measured, and it thus seems more difficult 
to obtain physical informations of the planet. Therefore, it 
is very important to extract them from the gravitational 
lensing effects. We have claimed that the rising spike, the 
sharp rise just before the ingress (or after the egress) of 
the transit, could be an excellent probe for finding planet 
parameters. It has a peculiar feature that the height has 
different dependence on the mass Ml, the radius Rl, and 
the orbital radius a from other observables, such as the depth 
and the duration of the transit. 

The rising spike can be large enough to be measured in 
high precision observation such as the Kepler, typically when 
the orbital radius is a few hundred au. It can be observed 
even if we take into account the limb darkening effect in 
spite of the fact that other peculiar features of the lensing 
effect on the transit light curve will be lost. 

Since the orbital radius is relatively large, one has at 
most one chance to observe in one's lifetime. The observ- 
ability of this kind of the transit could be estimated opti- 
mistically as 

7-, ^obs 



Torbit O, 



4 yr 



X 10"" ~ 3 X 10" 



(22) 



3 X 103 yr 200 au 

where the first factor is the chance to meet the transit, the 
second is the probability of the edge-on to the line of the 
sight, and the last is the number of the stars which will be 
searched in the Kepler in four years of its operation. Al- 
though it seems unlikely to find one by Kepler, one may 
observe in the future by using deeper and larger survey. 
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APPENDIX A: THREE CASES OF THE 
MICROLENSING AND THE OCCULTATION 

Here we complete the cases of the microlensing and the oc- 
cultation by the lens planet when the centres of the planet 
and the star are aligned. The total magnification is given by 
Eq.([6}, and its denominator is obtained as Eq.((7)|. There are 
three cases for the microlensing and the occultation: (I) Full, 
(II) partial, and (III) no occultations, as shown in Fig. \M 
Full occultation occurs when 9l > Oe, where both + and 
— paths are obstructed, while only — path is shielded in 
the outer region. Therefore, the numerator of Eq.® can be 
written as 



A+{l3)l3dl3d(j) 



TV 

2 



el-i3l + Os^el + 49| - Pl^JpI + 40 



,(A1) 



where Pl = P{Ol)- In the case (II), only partial occultation 
takes place in the outer region, where I3[—6l) < P < 0s, so 
that the numerator of Eq.© is given by 



27r /./3{-fli) 



A{l3)l3df3d(l) + 



A+{l3)l3dl3d(f> 



el~i3l + es Vol + 4e| - i3l + ^^l 



.(A2) 



In the last case, there is no occultation, which is just the 
same situation as the point mass lens. Thus, the numerator 
of Eq.(|6]) becomes 



A{l3)l3dl3d(f> = -KOs^Jel + m 



(A3) 



Therefore, the magnification is obtained as 




Figure Al. Schematic view of the microlensing and the occul- 
tation for three cases: Full, partial, and no occultations from the 
top to the bottom panels, respectively. 



in Fig. IBll while A and r in Fig. IB2|) have similar depen- 
dence on the physical parameters. Notice that one of these 
parameters (A in Fig. IBll while r in Fig. IB2[I changes very 
rapidly with respect to the parameter perpendicular to the 
two shown in the figures. 



1 



/3! + 



Am = 



(A4) 



(A5) 



where the tilde denotes the variable normalised with respect 
to Os such that /3l = Pl/Os and Oe = Oe/Os- 



APPENDIX B: OTHER SECTIONS OF THE 
ISOSURFACES OF THE RISING SPIKE, THE 
DEPTH, AND THE DURATION 

The rising spike ^, the depth A, and the duration r of the 
transit are functions of three variables: the mass Ml, the 
radius Rl, and the orbital radius a of the planet. It might 
be better to draw these isosurfaces in three dimensions, but 
may be difficult to look at. Therefore, we show some con- 
stant Ml sections in Fig. [5] because they seem to have very 
different curves in those sections. Here we show constant 
Rl and a sections for completeness in Figs. iBll and IB2l re- 
spectively. In these figures, two of the observables (^ and A 
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Figure Bl. Sections of tlie isosurfaces of ttie rising spike § 
(red), tlie depth A (green), and the duration r in days (blue) 
for Rl/Rj = 3, 1, and 0.3, from the top to the bottle panels, re- 
spectively. Shaded region represents where the hollow disappears 
A < 0, which we do not consider. 
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Figure B2. Sections of the isosurfaces of the rising spilte ^ (red), 
the depth A (green), and the duration r in days (blue) for a = 
300, 200, 100, and 10 au. Shaded region represents where the 
hollow disappears A < 0, which we do not consider. 



